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Abstract 

The linear response of non-equilibrium systems with Markovian dynamics satisfies a 
generalized fluctuation-dissipation relation derived from time symmetry and antisymmetry 
properties of the fluctuations. The relation involves the sum of two correlation functions of 
the observable of interest: one with the entropy excess and the second with the excess of 
dynamical activity with respect to the unperturbed process, without recourse to anything 
but the dynamics of the system. We illustrate this approach in the experimental determi- 
nation of the linear response of the potential energy of a Brownian particle in a toroidal 
optical trap. The overdamped particle motion is effectively confined to a circle, undergoing 
a periodic potential and driven out of equilibrium by a non-conservative force. Independent 
direct and indirect measurements of the linear response around a non-equilibrium steady 
state are performed in this simple experimental system. The same ideas are applicable 
to the measurement of the response of more general non-equilibrium micron-sized systems 
immersed in Newtonian fluids either in stationary or non-stationary states and possibly 
including inertial degrees of freedom. 

1 Introduction 

The linear response of systems in thermodynamic equilibrium is generally described by the 
fluctuation-dissipation theorem [T|. It provides a simple relation between the equilibrium fluc- 
tuations of an observable Q with the response due to a small external perturbation h s changing 
the potential at time s as U — > U — h s V: 

RQv(t-s)=P-^-(Q(t)V(s)) . (1) 

In equation (TjQ) B,Qv(t — s) = S(Q(t))h/Sh s \h=o is the linear response function of Q at time t > s; 
(Q(t)V(s)}o is the two-time correlation function between Q and V measured at equilibrium; 
the brackets (. . .)h denote the ensemble average in the state perturbed by h s so that (. . .) 
corresponds to the ensemble average at equilibrium (h s = 0). The inverse temperature of 
the equilibrium system, f3 = l/fc^T, appears as a multiplicative factor. Hence, equation ((T|) 
represents a useful tool in experiments and simulations to explore indirectly the linear response 
regime from fluctuation measurement completely performed at thermal equilibrium. Vice versa, 
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one can obtain information on microscopic fluctuations from non-equilibrium measurements of 
response functions or susceptibilities by applying sufficiently weak external fields. 

In general, equation ([T]) fails to describe the linear response of systems already prepared 
in a non-equilibrium state. This situation is relevant in real mesoscopic systems that usually 
operate far from equilibrium due to cither non-conservative/time-dependent forces exerted by 
the experimental apparatus or external flows and gradients applied at the boundaries. For 
instance, the developcmcnt of micro and nano techniques {e.g. optical tweezers and atomic 
force microscopes) has allowed one to mechanically manipulate colloidal particles, living cells 
and single molecules of biological interest with forces ranging from pN to fN. In this kind 
of experiments, thermal fluctuations can be comparable or larger than the typical external 
perturbations necessary to determine Rqv ■ Then, for these systems it is more reliable in practice 
to measure non-equilibrium fluctuations than linear response functions. 

On the theoretical side, several works have recently dealt with the problem of the exten- 
sion of the fluctuation-dissipation theorem around non-equilibrium steady states [H OH El 03 G3 
[7J |H1 [HJ EH HI]- In most of the formulations, an additive extra term on the right-hand side 
of equation ([T]) appears as a non-equilibrium correction of the fluctuation-dissipation relation 
due to the broken detailed balance. Different physical interpretations of the corrective term 
are provided in the literature. We specially highlight the roles of the non-vanishing probability 
current [5] and the conjugate variables to the total entropy production [SJ [TT] in the analyti- 
cal expression of the corrective term, leading independently to a Lagrangian interpretation for 
Langevin systems with first-order Markovian dynamics. Following these two approaches an 
equilibrium-like fluctuation-dissipation relation can be restored in the Lagrangian frame of the 
local mean velocity of the system for the right choice of observables [6) [12] . An alternative 
interpretation, the entropic-frenetic approach, has been proposed in terms of time-symmetric 
and time-antisymmetric properties of the fluctuations |13 1 I14 [ [T5|. Unlike previous formulations, 
the entropic-frenetic approach does not involve the stationary probability distribution of the 
degree of freedom of interest but only explicit observables. In addition, it holds even in more 
general situations including non-steady states and inertial degrees of freedom provided that the 
dynamics is Markovian |15j . Therefore it is more accessible from the experimental point of view. 

The present paper provides the first experimental investigation of the modified fluctuation- 
response relation of [13111 4 j in a simple experimental system driven into a non-equilibrium steady 
state: a Brownian particle whose overdamped motion is effectively confined to a circle by a 
toroidal optical trap [TB] . As the new formula relating response and fluctuations is exact and does 
not involve further notions such as effective temperatures, it can be tested at full value. Moreover, 
we do not need the explicit knowledge (analytic form) of the stationary distribution. Hence 
the experimental linear response of the system can be safely determined by two independent 
methods. In the first method one control parameter of the non-equilibrium steady state is 
physically perturbed and the corresponding response of the particle is directly measured. In the 
second method, only the unperturbed non-equilibrium fluctuations of the position of the particle 
are measured and the same linear response is determined based on the generalized fluctuation- 
dissipation formula and on a suitable model of the dynamics. All of this is done here under 
very well controlled experimental conditions, so that all can be checked explicitly and every step 
can be analyzed separately. That control allows us to understand the usefulness of the new 
formalism so that it will be able to be employed also in natural situations where such a control 
is really absent, and e.g. the driving or potential are not exactly known. 

The paper is organized as follows: in section [2] we briefly present the generalized approach 
to linear response based on the entropic-frenetic formulation of a fluctuation-dissipation rela- 
tion for systems with Markovian dynamics and the main goals of our experimental study with 
respect to this formalism. In section [3] we describe the main features of the experiment and the 
ID Langevin model of the translational motion of the particle in the toroidal trap. In section [4] 



2 



we present the results of the two methods to measure the linear response function of the sys- 
tem. We show that within our experimental accuracy both methods lead to the same values of 
the linear response function when taking into account the corrections given by the generalized 
fluctuation-dissipation relation. Then, we discuss their technical limitations and advantages 
from the experimental point of view. We also depict a simple example in order to show the 
flexibility of the generalized fluctuation-dissipation formula: once the unperturbed fluctuations 
of the proper degree of freedom are measured, one can readily compute the linear response of the 
system under more complex time-dependent perturbations. Finally we present the conclusion. 



2 Generalized approach to linear response 

We briefly present the formulation of the fluctuation-dissipation relation developed in p~3l [T4l [T5] 
for the special case of stochastic systems described by a finite number of degrees of freedom {q} 
with overdamped Markovian Langevin dynamics, in presence of a potential U(q). We consider 
a system in contact with a thermal bath at temperature T and driven into a non-equilibrium 
steady state by a non-conservative force. We focus on the average value of an observable Q(q) 
at time t, denoted by (Q{qt))o- We are also interested in the mean value (Q(qt))h of Q(q), when 
a small time dependent perturbation h s is applied to U(q), namely U(q) — > U(q) — h s V(q). As 
formally shown in [T21H1], (Q(qt))h is given at linear order in h s by 

(Q(qt)) h = (Q(qt)h + I Rgv(t, s)h s ds, (2) 

J ~ oo 

where the linear response function Rqy obeys the generalized fluctuation-dissipation relation 

= ~ s (V(q s )Q(q t )) - ^(LV(q 8 )Q(q t )) . (3) 

In equation ((3j) , L is the generator of the unperturbed Langevin dynamics which determines the 
time evolution of any single-time observable O(q): d(0(q t ))o/dt = ((LO)(q t )) . 

It should be noted that RQv(t,s) is operationally obtained by applying an instantaneous 
delta perturbation at time s and measuring (Q(qt))h — (Q(qt))o- I n experiments it is always 
more reliable to apply a Heaviside perturbation (h s = for s < 0, h s = h = const, for s > 0) 
instead. This procedure directly yields the integrated response function 

,.s [* D ( . x A {Q(qt))h - {Q(gt))o 

XQv{t)=J R Q v{t,s)ds = , (4) 

defined over the time interval [0,i]. Therefore, in the following we consider the integral form of 
equation ^ 

XQv(t) = £[C(t)+K(t)] t (5) 

where the term 

C(t) = (V(q t )Q(q t ))o - (V(q )Q(qt))o, (6) 

can be interpreted as a correlation between the observable Q(qt) and the excess in entropy 
produced by the Heaviside perturbation during the interval [0,t]: [hV(qt) — hV(qo)]/T. On the 
other hand, the term 

K(t) = - [ (LV(q s )Q(q t )) ds, (7) 
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can be interpreted as minus the correlation between Q{qt) and the integrated excess in dynamical 
activity or frenesy: (3 f Q LV(q s )hds, which quantifies how frenetic the motion is due to the 
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perturbation with respect to the unperturbed process. The frenesy (3LV(q) can be regarded as a 
generalized escape rate of a trajectory from a given phase-space point q. In references [13lfT4lfT5] 
it has been widely discussed that the origin of the entropic C(t) and the frenetic K{t) terms can 
be traced back to time-antisymmetric and symmetric properties of the fluctuations, respectively. 

The non-equilibrium fluctuation-response formulae @ and ([3]) result directly from the theory 
of the Langevin equation, and hence, the validity of a Langevin model for our experiment would 
thus also imply this particular fluctuation-response formula. Since for the present experiment 
the Langevin equation has been largely verified before, thus implying in theory the response ([2])- 
(0, the point of the present experimental study is not in the very first place to validate 
although we still want to make sure it is exact. Rather, the point of the present investigation is 
much more twofold: 

1) to show how the (very new) frenetic term ([7]) can indeed be measured independently and 
directly, without recourse to anything but the dynamics. In particular, no knowledge of 
the analytical expression of the stationary distribution is used; 

2) to show how the new formalism ([2])-([7|) of response can be used to determine unknown 
forces, both conservative and non-conservative. For example, when the formula is accepted 
as such, the correction to the equilibrium formula gives useful information about how the 
breaking of detailed balance gets established. Alternatively, by changing the driving we 
can discover the potential that acts on the particle, etc. 

These points are addressed in the following sections under proper experimental control and data 
analysis, which allows us to understand how the generalized fluctuation-response formalism can 
be applied to a real non-equilibrium experiment. 

3 Colloidal particle in a toroidal optical trap 

3.1 Experimental description 

We specifically study the linear response of a single colloidal particle driven out of equilibrium 
by a non-conservative force in presence of a non-linear potential. We recall the main features 
of the experiment, previously described in detail in reference |16j . where the same experimental 
set-up was used in the context of a modified fluctuation-dissipation relation with a different 
interpretation from the one described in the present paper. In our experiment the Brownian 
motion of a spherical silica particle (radius r = 1 fim) immersed in water is confined on a 
thin torus of major radius a = 4.12 fim by a tightly focused laser beam rotating at 200 Hz (see 
figure HJa)). The rotation frequency of the laser is so high that it is not able to trap continuously 
the particle in the focus because the viscous drag force of the surrounding water quickly exceeds 
the optical trapping force. Consequently, at each rotation the beam only kicks the particle a 
small distance along the circle of radius a. During the absence of the beam (w 5 ms), the 
particle undergoes free diffusion of less than 40 nm in the radial and perpendicular direction to 
the circle. Thus, the particle motion is effectively confined on a circle: the angular position 9 
of its barycenter is the only relevant degree of freedom of the dynamics. In addition, a static 
light intensity profile is created along the circle by sinusoidally modulating in time the laser 
power which has a mean value of 30 mW and a modulation amplitude of 4.2 mWpp at the 
same frequency as the rotation frequency of the beam. Figure QJa) sketches this experimental 
configuration on a snapshot of the colloidal particle in the toroidal trap. The water reservoir 
acts as a thermal bath at fixed temperature (T = 20 ± 0.5°C) providing thermal fluctuations to 
the particle. The viscous drag coefficient at this temperature is 7 = 1.89 x 10~ 8 kg s . 
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Figure 1: (a) Snapshot of the colloidal particle in the toroidal optical trap. The vertical arrow 
indicates the position 9 = whereas the curled arrow shows the direction of the rotation of the 
laser beam, (b) Experimental potential profile (solid line) and probability density function of 
9 (dashed line) for the non-equilibrium steady state generated by the non-conservative force /. 
Inset: Typical steady-state trajectory {9t, < t < 66.67 s} used to compute p(9). 



3.2 Model 

For the experimentally accessible length and time scales the dynamics of 9 is modeled by the 
first-order Langevin equation 

B = -Ati/(0)+F + t, (8) 

as extensively verified in [TBI [T71 HB HH HD] ■ A(f)(0) is a periodic non- linear potential [A(f>(9) = 
Acf)(9 + 27r)] of amplitude A = 0.87 rad 2 s _1 created by the periodic modulation of the laser 
intensity. The normalized angular profile (f)(6) of the potential (max{| </>(#) |} = 1) is a slightly 
distorted sine function as a result of unavoidable experimental static defects of the toroidal 
optical trap (e.g. optical aberration). F = 0.85 rad s _1 is a constant force acting in the direction 
of the laser rotation which is associated to the mean kick of the beam. £ is a white noise process of 
zero mean and covariance = 2DS(t~s) with bare diffusivity D = k B T/(~fa 2 ) = 1.26 x 10~ 2 

rad 2 s _1 , which models the thermal fluctuating force exerted by the water molecules. F is non- 
conservative (Jq™ Fd9 = 2nF > 0) since the motion takes place on a circle, driving the system 
out of equilibrium. The physical non-conservative force and the potential are / = jaF = 66 
fN and U(9) = ja 2 A(j>(0) = 68. 8k bT4>(0), respectively. The experimental potential profile 
U(9) is plotted as a continuous black line in figure (Hb). Note that at thermal equilibrium 
(F = 0) the particle motion would be tightly confined around the potential minimum with the 
stochastic variable 9 distributed according to the Boltzmann density p e q(9) oc exp[— /3U(6)]. 
However, due to the thermal fluctuations and the non-conservative force F > 0, the particle 
is able to go beyond the potential barrier and explore the whole circle. In the non-equilibrium 
situation with constant F, A,D > 0, the angular position 9 settles in a stationary probability 
density p(9) ^ Peq(@) that admits an analytical expression found in [21] , The experimental 
non-equilibrium density p(9) is shown in figure HKb) as a dashed line. A constant probability 
current j = (0(t)) o /(2w) = [F — A<p'(6)]p(0) - Dd e p(9)\ > 0, in the direction of F appears 
reflecting the broken detailed balance of the dynamics. For the experimental conditions one 
finds j = 3.76 x 10~ 2 s _1 corresponding to a mean rotation period of 26.6 s for the particle. 
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Both p(9) and j, determined from 200 independent experimental time series of the angular 
position of the particle {9 t , < t < 66.67 s}, allow one to precisely compute the values of / 
and U(9), as described in detail in |T51 |T7] . One example of such a steady-state time series is 
depicted in the inset of figure [TJb). 

4 Measurement of the linear response function around a 
non-equilibrium steady state 

Now we proceed to determine the linear response of the particle motion when slightly perturbing 
the non-equilibrium steady state previously described. For experimental simplicity we consider 
a step perturbation to the potential amplitude A — > A + 6 A, so that the perturbation and its 
conjugate variable are h = —5A and V(q) = 4>{9), respectively. It should be noted that the 
dynamics of 9 is strongly nonlinear as the particle undergoes the periodic potential A<p(9) but 
for sufficiently small values of 5 A the response of <f>(9) can still be linear. In reference [T5] , 
where we used the same experimental data as in the present paper, we have already checked 
that for a perturbation SA of the potential amplitude A in the range |<L4| < 0.07A the linear 
response regime holds. This was done by computing the experimental integrated response x(t) 
of the observable Q(9) = sin 9 for two different values of SA: 0.05A and 0.07A We showed 
that within the experimental error bars x(t) is the same in both cases. The independence of 
X(i) with respect to 5 A experimentally demonstrates that the system is actually in the linear 
response regime at least for |<L4| < 0.07A and times < t < 3.3 s after the application of the 
Heaviside perturbation. 

In order to illustrate the meaning of the non-equilibrium linear response relation of equation 
in this case, we perform two different kinds of independent measurements: direct and indirect, 
as explained in the following. 



4.1 Direct measurement of the linear response function 

First, we consider the direct measurement of the integrated response function xqv for the 
periodic observable Q(q) = 4>(9) = 4>(9 + 2ir). This observable times 7a 2 A represents the 
instantaneous potential energy of the particle. This is experimentally accomplished by applying 
a step perturbation to the amplitude of the sinusoidal laser power modulation, 4.2 mWpp — > 
4.4 mWpp, but keeping constant the power offset at 30 mW so that F remains constant. This 
results in a perturbation of the potential amplitude 6 A = 0.05A. Then, the integrated response 
function of (f)(9) at time t > due to the perturbation — 5 A applied at time is given by 

XQv{t) = , (9) 

where {...)sa an d (•••)() denote the ensemble averages of the perturbed 9 t ^sA and unperturbed 9t 
trajectories, respectively. In order to decrease the statistical errors in comparison of the terms 
in equation ([9]), for a given perturbed trajectory 9t : sA we look for as many unperturbed ones 
9t as possible starting at time t* such that 4>{9f) = 0(#o.<5a)- Then we redefine t* as t = in 
equation ([9]), as depicted in figure [2ja). The unperturbed trajectories found in this way allow 
us to define a subensemble over which the average (...)o is computed. On the other hand, the 
average {...}sa is performed over 500 independent realizations of SA that are enough for a fair 
determination of the integrated response. The resulting curve XQv(t) and its corresponding 
error bars ±<r x (i) as functions of the integration time t starting at the instant of the application 
of the step perturbation are shown in figures EJc) and EJd) . 
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Figure 2: (a) Examples of a perturbed trajectory (red thick dashed line) and three unperturbed 
steady-state trajectories (thin solid lines) used to compute the integrated response function given 
by equation Inset: expanded view at short time. The unperturbed trajectories are chosen to 
start at the same value as the perturbed one. (b) Potential profile <fi(9) locally fitted as a third- 
order polynomial (pLocifi) around each value of 9. Their derivatives are computed from 4>Loc(9)- 
(c) Integrated linear response function XQv(t), entropic C(t) and frenetic K(t) terms and the 
corresponding indirect measurement of the response [C(t) + K(t)]/2 for Q = <p(0), h = —SA and 
V = 4>{0), as functions of the integration time t. (d) Expanded view of (c). 
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4.2 Indirect measurement of the linear response function 

The same response information can be obtained indirectly from correlation measurements of 
the unperturbed non-equilibrium steady-state fluctuations (8 A = 0) of 9 t when properly using 
equation ([5]) . For the Langevin dynamics of 9 described by equation the analytical expression 
of the generator L is 

L = (F-Atj/(0))do + D8%. (10) 
Hence, in this case equation ([5]) reads 

DXQv(t) . <m±m, 

where the entropic and frenetic terms are 

C(t) = {4>(p t )<K9 t ))o - (^)^)) , (12) 

K(t) = - [ ds{[D<j>"{e s ) + {F-A4>'{e s ))4>'{e s )\<p{9 t )) Q , (13) 

Jo 

respectively. At this point it is clear that we need to know the potential profile 4>{0) for the 
indirect method. The integrand of equation (|13[) involves the instantaneous values of <f>(9) and 
its derivatives 4>'{6) and 4>"(9). In order to take into account the non-sinusoidal distortion of the 
potential profile, we perform a local polynomial fit <f>Loc of cj> around each value of 9 £ [0,27r). 
Then the instantaneous value of the observable 4>(9 t ) at time t is approximated by 4>Loc(9t) 
either for an unperturbed or a perturbed trajectory. The local polynomial approximation (j>Loc 
and its derivatives </>^ oc , fiLoc arG plotted in figure |5Jb) showing the non-sinusoidal distortion. 

The resulting curves C(t) and K(t) as functions of the integration time t are plotted in figure 
[2{c). At thermal equilibrium (F = 0) one should find that C(t) = K(t) for all t > because 
of the time reversibility and stationarity of the two-time correlations leading to the equilibrium 
fluctuation-dissipation relation D\Qv(t) = C{t). On the other hand, in the present case K(t) 
reaches negative values of the same order of magnitude as the positive values of C(t). This 
reflects the experimental conditions far from thermal equilibrium of the system. The curve for 
K(t) represents the first experimental result concerning the direct measurement of the dynamical 
activity along a trajectory [13l HU [15] . The average of these two quantities [C (t) + K (t)] / 2 is one 
order of magnitude smaller. This average, which is an indirect measurement of the integrated 
response function according to equation (|11|). agrees very well with the direct measurement 
of xqv within the experimental error bars, as shown in figure [U(d). All is consistent with 
the results of [16], see [6j [13], but the experimental approach here is quite different. The 
experimental verification of previous equivalent generalized fluctuation-dissipation relations for 
the same experimental system as the one studied here [T6l [18, 20J involves the local mean velocity 
of the particle (i.e. the ratio between the steady-state current and the probability density). Here 
we measure directly explicit correlation functions C(t) and K(t) without recourse to and indeed 
without need for the expression for the stationary distribution. Note that the experimental 
density profile p(9) (and not its analytical expression) is used in practice only for calibration 
purposes [HI [17] to determine the experimental parameters A and F. However, once these 
parameters are known the data analysis necessary to compute C(t) and K(t) completely relies 
on the dynamics, i.e. the measurement of the time series 9f . Thus, this procedure is suitable 
for the study of the linear response of more complex Markovian non-equilibrium situations (e.g. 
multiple degrees of freedom, non-stationary states, several thermal baths) where local mean 
velocities are not easily measurable in experiments. 
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4.3 Discussion 



We remark that for this kind of micron-sized system the relation ([5]) actually represents a feasible 
indirect method to access the linear response regime far from thermal equilibrium. This is 
because all the parameters of the unperturbed dynamics are known a priori or can be determined 
in situ without any external perturbation of the non-equilibrium state. On the other hand, the 
direct measurement of the linear response function exhibits a number of technical difficulties in 
practice. First, a vanishingly small Hcavisidc perturbation —hV(q) to the initially unperturbed 
potential U (q) is ideally required. Otherwise spurious effects quickly bias the measurement of 
XQV, specially when the system is strongly non- linear. Second, one requires an extremely large 
number of independent realizations of h to resolve XQV as the perturbation — hV(q) must be 
chosen very weak, typically smaller than the thermal fluctuations of the energy injected by the 
environment. This is evident on the results of the direct measurement of XQV for the colloidal 
particle, see figure [^d). For integration times t < 3 s the agreement between DxQv(t) an d 
[C(t) + K(t)}/2 is excellent. Then, for t > 3 s the finite sampling leads to deviations between 
the two methods and increasingly large error bars for the direct measurement of XQV- These 
drawbacks are skipped when implementing the indirect method measuring the unperturbed 
quantity [C(t) + K(t)]/2. In addition, for a steady state like the one experimentally studied 
here one can improve dramatically the statistics by performing an additional time average over 
a window [0, t max ] : C(t) = J* mam C(t + u)du/t max , K(t) = J mam K(t + u)du/t max , as actually 
done for the curves C, K and (C + K)/2 in figurcs^c) and^d). However, one must be careful 
when performing the time average. This is because the value of t max may affect the resulting 
values of [C(t) + K(t)]/2 as t increases, specially for correlation functions involving strongly 
fluctuating quantities such as velocities, as recently discussed in [20] . For the curves shown in 
figures GJc) and[2jd) we verified that their shapes are not significantly influenced by t max . 

4.4 Application example 

Finally, we illustrate the usefulness of the indirect measurement of the linear response function 
to study the temporal behavior of the mean potential energy of the particle (U (&t))h under small 
external perturbations h s more intricate than a simple Hcavisidc function. This is done with the 
same unperturbed experimental data used in subsection 14.21 without carrying out the different 
physical realizations of h s . We concentrate on a sinusoidal perturbation starting at time s = 

h s = h sin27r/ s. (14) 

either to the phase or to the amplitude of the potential around the steady state. First, we consider 
the case of a small phase perturbation a s = ao sin27r/os with ag <C 1 : A<fi(6) — > A<p{9 + a s ) « 
A(j){ff) + A(j)'{ff)a s so that ho = — Aao in equation (|T4|) and V = 4>'{0). One must compute first 
the integrated response function of Q = cj)(9) given by [C{t) + K(t)]/(2D) by inserting the right 
V and Q in equations © and 0. The resulting curves C(t), K(t) and \C{t) + K(t)]/(2D) 
are plotted in figure EJa) . Next, using equation ([2]) the experimental impulse response function 
[dtC(t — s) + d t K{t — s)]/(2D) must be convolved with h s given by equation (fT4")) . In this way 
one finds that the mean potential energy of the particle oscillates around the non-equilibrium 
steady state value (U(9))o = —4.71cbT as shown by the dashed blue line in figure [3(b) for 
ao = 0.05 rad and fo = l Hz. The oscillations exhibit a delay time (At 0.23 s) with respect to 
a s and a slow transient (~ 15 s) corresponding to the decay of the non-equilibrium stationary 
correlations. As t increases the oscillations settle around (U(9))q with a constant amplitude 
AU ~ Q.2ksT. Now, we consider a sinusoidal time-dependent perturbation to the potential 

lr The only restriction in the choice of t max is that the ensemble of averaging intervals [So,St ma:c ] must cover 
the whole circle to sample correctly the steady state. 
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Figure 3: a) Integrated response function of the observable Q = (f)(9) for a small perturbation 
of the potential phase as a function of the integration time t. Inset: expanded view of [C(t) + 
K(t)]/2. (b) Sinusoidal timc-dependent perturbation —h s (solid black line) of the static potential 
A(f>. Resulting mean potential energy of the Brownian particle for a phase perturbation (dashed 
blue line) and an amplitude perturbation (dotted red line) for — ho/ A = 0.05 and fo = 1 Hz. (c) 
Asymptotic values of oscillation amplitude of the potential energy and (d) the delay time with 
respect to — h s for each kind of perturbation. The black dashed lines represent the values that 
would be obtained around thermal equilibrium (F = 0), given by equations (|17[) and (|18l) . 
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amplitude: SA S = 5Aq sin27r/os with the same strength (—ho/ A = 0.05) and frequency as 
before. In this case ho = — 5Ao and V = (j>{6). Following the same procedure with [C + K]/ (2D) 
shown in figure [2]Jd), one finds a different qualitative behavior of (U(9 t ))h, as depicted by the 
dotted red line in figure [21b) . At the beginning the mean potential energy responds in the 
opposite direction to SA S . Then, as t becomes larger than the slow non-equilibrium transient 
(U(9))h oscillates around (U(9)) with a constant amplitude AU « O.bksT and a delay time 
At rj 0.26 s. For both types of perturbations one can write the asymptotic dependence of 
(U(0 t ))h on t > 15 s as 

(U(6 t )) h = (U(6 t ))o ± A?7sin[27r/o(i - At)], (15) 

where the positive and negative signs stand for the phase and amplitude perturbations, respec- 
tively. The values of AU and At depend on the frequency fy. In figures [3Jc) and^d) we show 
this dependence. We now compare these far-from-equilibrium results with those that would be 
obtained when applying h s around thermal equilibrium (F = 0). In such a case the particle 
motion is tightly confined to the harmonic part of the potential around the minimum 9 m — 3tt/2: 
4>(9) m — 1 + (6 — 9 m ) 2 /2. After some algebra using this approximation one finds the expression 
for (U(9 t ))h when perturbing thermal equilibrium 

(U(9 t )) h = (U(6 t ))o±AU{sm[2irf (t-At)} + 

e - 2At sm27rf a At}, (16) 



where (U(9 t )) = -68.3fc B T and 



AU = -^ ^ (17) 

A 2(1 + tt 2 /oA4 2 ) 1/2 ' ' ' 

At = — !— arctan ( ^ ) , (18) 



2nf V A 

either for a phase (positive sign) or an amplitude (negative sign) perturbation. Note that for 
t 3> (2 A) -1 , equation (|T6"|) exhibits the same qualitative behavior as (fTS]). We plot the curves 
given by eqsuations ([TT1) and (jTSJ) in figures ^c) and|3Id), respectively, for the same values of 
the parameters ho and A as before. Unlike the behavior close to equilibrium, the oscillation 
amplitude AU strongly depends on the perturbed parameter around the non-equilibrum steady 
state: it is more sensitive to amplitude perturbations than to phase perturbations. In addition, 
the far-from-equilibrium values are two orders of magnitude larger than that given by equation 
(TIT)) . By contrast, the delay time At is not significantly affected by the far-from-equilibrium 
nature of the system. It is almost independent of F and of the type of perturbation and it 
converges to equation (fTS)) as fo increases. 



5 Concluding remarks 

We have experimentally studied the non-equilibrium linear response of the potential energy 
of a Brownian particle in a toroidal optical trap. We gain insight into the application of the 
fluctuation-dissipation relations far from thermal equilibrium in this non-linear system with a 
single relevant degree of freedom. In particular, we show that the entropic-frenetic approach 
is appropriate and feasible for the study of the linear response of non-equilibrium micron-sized 
systems with a small number of degrees of freedom immersed in simple fluids. Non-trivial 
linear response information can be obtained from purely unperturbed measurements of the non- 
equilibrium fluctuations provided that the parameters describing the dynamics are known. Our 
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experiment reveals that the indirect determination of the linear response function is less time- 
consuming, more accurate and more flexible than the direct perturbation of the non-equilibrium 
system. Similar ideas are expected to be applicable to more complex micron-sized systems 
such as atomic-force microscopy experiments and ensembles of colloidal particles in simple non- 
equilibrium conditions where local mean velocities are difficult to measure. 
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